A relativistic type Ibc supernova without a detected 
7— ray burst 
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Long duration gamma-ray bursts (GRBs) mark^ the explosive death of some 
massive stars and are a rare sub-class of Type Ibc supernovae (SNe Ibc). They 
are distinguished by the production of an energetic and collimated relativistic 
outflow powered^ by a central engine (an accreting black hole or neutron star). 
Observationally, this outflow is manifested'^ in the pulse of gamma-rays and a 
long-lived radio afterglow. To date, central engine-driven SNe have been dis- 
covered exclusively through their gamma-ray emission, yet it is expected^ that 
a larger population goes undetected due to limited satellite sensitivity or beam- 
ing of the collimated emission away from our line-of-sight. In this framework, 
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the recovery of undetected GRBs may be possible through radio searches^'^ for 
SNe Ibc with relativistic outflows. Here we report the discovery of luminous ra- 
dio emission from the seemingly ordiuciry Type Ibc SN 2009bb, which requires 
a substantial relativistic outflow powered by a central engine. The lack of a co- 
incident GRB makes SN 2009bb the flrst engine-driven SN discovered without 
a detected gamma-ray signal. A comparison with our extensive radio survey of 
SNe Ibc reveals that the fraction hcirboring central engines is low, ~ 1%, mea- 
sured independently from, but consistent with, the inferred^^ rate of nesirby 
GRBs. Our study demonstrates that upcoming optical and radio surveys will 
soon rival gamma-ray satellites in pinpointing the nearest engine-driven SNe. 
A similcir result for a different supernova is reported^ independently. 

On 2009 Mar 21.1 UT, the Chilean Automatic Supernova Search Program (CHASE; 
Ref. 31) discovcrcd^'^ a bright optical transient through repeated imaging of the nearby 
spiral galaxy NGC 3278 at a distance, c? 40 Mpc. The new object was offset 22 arcsec 
(4.2 kpc) from the center of the galaxy and located within its star-forming disk. Optical 
spectroscopy obtained on Mar 28.1 UT revealed^^ that the transient was a young Type Ibc 
SN (hereafter SN2009bb) lacking evidence for Hydrogen in the explosion debris. Based 
on the previous non-detection of SN 2009bb on Mar 19.2 UT, we tightly constrain the SN 
explosion date to be Mar 19 ± 1 UT (see Suppl. Info.). 

Using the Very Large Array (VLA) on Apr 5.2 UT, we discovered a coincident radio 
counterpart at a(J2000) = 10'^3r^33.87'^ and (5(J2000)=-39°57'30.1" (±0.7 arcsec in each 
coordinate) and with a flux density, F^, — 24.53 ± 0.06 mJy, at frequency, u — 8.46 
GHz. This corresponds to a spectral radio luminosity oi L,^ ^ 5 x 10^^ erg Hz~^ at 
At 17 days after explosion, more luminous than any other SN Ibc observed^'^'^^ on 
a comparable timescale. Instead, the radio properties of SN2009bb arc consistent with 
the sample of nearby {z < 0.1) GRBs, observed to consistently yield^^ lower relativistic 
energies than "classic" GRBs preferentially discovered at larger distances. Further VLA 
observations of SN2009bb revealed a power-law flux decay, F^^sag ghz ~ t~^ '^, in line with 
the radio afterglow evolution seen^^ for the nearest gamma-ray burst, GRB 980425 at a 
similar distance of o? « 38 Mpc (Figure 1). 
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Unlike the optical emission from SNe which traces only the slowest explosion debris, 
radio observations uniquely probe^^ the fastest ejecta as the expanding blastwave (velocity, 
v) shocks and accelerates electrons in amphfied magnetic fields. The resulting synchrotron 
emission is suppressed by self- absorption (SSA) producing a low frequency radio turn- 
over that defines the spectral peak frequency, Up. Combining our observations from the 
VLA and the Giant Meterwave Radio Telescope (GMRT), the radio spectra of SN2009bb 
are well described by an SSA model across multiple epochs (Figure 2). Prom our earliest 
spectrum on Apr 8 UT {At ^ 20 days), we infer Up ^ 6 GHz and a spectral peak luminosity, 
L^,p ^ 3.6 X 10^8 erg s-^ Rz'K 

Making the conservative assumption that the energy of the radio emitting material 
is partitioned equally into accelerating electrons and amplifying magnetic fields (equipar- 
tition), the properties of the SSA radio spectrum enable^^'^^ a robust estimate of the 
blastwave radius, i? 2.9 x W^^{L^^p/10'^^ erg s"^ liz-^f/^^{up/5 GHz)-i cm. Luminous 
synchrotron sources with a low spectral peak frequency thus require larger sizes (Figure 3). 
For SN 2009bb, we infer R rs 4.4 x 10^^ cm at At ~ 20 days and thus the mean expansion 
velocity is R/At — 0.85 ± 0.02c, where c is the speed of fight. The transverse expansion 
speed, T/Sc = R/At indicates that the blastwave is relativistic, F > 1.3, at this time [bulk 
Lorentz factor F = (1 - /32)-i/2 ^j^h /3 = y/c]. This is a lower limit on the initial velocity 
since the radio evolution indicates that the blastwave decelerated early on. We further find 
that the radio emission requires a minimum energy, E — (1.3±0.1) x 10^^ erg, coupled 
to the relativistic outfiow and comparable to the values inferred^^'^^'^^'^^ from the radio 
afterglows of nearby GRBs (see Figure 4). 

These conclusions are robust; the blastwave velocity is insensitive^^'^^ to deviations from 
equipartition while the relativistic energy can only be higher. In view of these constraints, 
we note that shock-acceleration in some SNe Ibc may^^ couple a minute fraction (< 0.01%) 
of the total energy, £^tot; to material with a trans-relativistic velocity. However, this scenario 
would require an exceedingly high total energy for SN2009bb, £Jtot ^ 10^^ erg, a factor of 
10^ higher than the total explosion energies inferred for SNe Ibc. We conclude that the 
energetic and relativistic outfiow from SN 2009bb was powered by another energy reservoir, 
a central engine. To date, engine-driven SNe have been discovered exclusively through their 
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gamma-ray emission, making SN2009bb the first to be identified by its long- wavelength 
signal. 

Motivated by our discovery of an engine-driven relativistic outflow, we searched for a 
gamma-ray counterpart in temporal and spatial coincidence with SN 2009bb. During our 
bracketed explosion date estimate, the all-sky Interplanetary Network (IPN; Ref. 44) of 
high energy satellites did not detect a coincident GRB (see Suppl. Info.). Based on the IPN 
sensitivity and detection efficiency, we place an upper limit on the gamma-ray fluence of 

< 5 X 10~^ erg cm""^, corresponding to an energy of Ej < 10^^ erg (band, 25-150 keV) if 
the relativistic outflow was isotropic. This hmit is a factor of two higher than the isotropic- 
equivalent observed''^ from GRB 980425, and thus it is possible that SN2009bb gave 
rise to a similar (albeit undetected) signal. At the same time, these limits cannot exclude 
scenarios in which the SN (i) powered a GRB directed away from our hne-of-sight, or (ii) 
did not produce any gamma-rays. SN 2009bb observationally demonstrates the limitation 
of using gamma-ray satellites as a primary tool to identify nearby engine-driven explosions. 

In this context, we note that our VLA observations of SN 2009bb were obtained as part 
of an extensive radio survey of 143 optically-discovered local SNe Ibc designed to recover 
relativistic SNe without detected gamma-ray counterparts. This systematic radio study 
revealed^^'^'^'^^ no additional relativistic SNe, instead indicating typical maximum blast- 
wave velocities for SNe Ibc of /3 ~ 0.1 (Figures 3 and 4). From this sample, we therefore 
estimate the fraction of engine-driven SNe to be just 0.7^.1'^% (Ic). This is measured inde- 
pendently from, and yet consistent with, the fraction inferred from the relative rates^^'^^'^*^ 
of nearby GRBs and all SNe Ibc, 1.41q;3% (Icr). Our long-term study thus confirms that 
engine-driven SNe are uncommon. 

The infrequency of relativistic outflows among massive star explosions implies that 
their progenitor stars share an essential and rare physical property. Observations of the 
explosion environment may offer unique clues. On a local scale (< 1 pc), where the 
environment was shaped directly by the evolution of SN 2009bb's progenitor star, we flnd 
evidence for a pre-explosion mass loss rate of M = (2.0 ± 0.2) x 10~^ yr~^, consistent 
with the wide distribution of values inferred^^'^^'^^ for nearby GRBs. However, the large- 
scale (> 1 kpc) environmental properties differ*^° from those of nearby GRB host galaxies. 
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showing evidence for a super-solar metallicity that exceeds the proposed^^'^^ cut-off for 
relativistic explosions. We conclude that there is a broad diversity in the environments 
of engine-driven explosions, and therefore host galaxy properties alone cannot be used to 
discriminate between ordinary and engine-driven SNe. 

With the advent ^''''^^ of wide-field optical surveys (e.g., Palomar Transient Factory, Pan- 
STARRS) the discovery rate of young, local SNe Ibc will effectively quadruple over the 
next 3-5 years (see Suppl. Info). Coupled with the ten-fold increase in the sensitivity of the 
Expanded Very Large Array (expected 2010; Ref. 51), relativistic SNe will be uncovered 
at an increased rate of ~ 1 per year within d < 200 Mpc. This is ~ 3 times higher than 
the rate at which nearby GRBs are discovered with current gamma-ray satellites. Thus, 
while such explosions have historically been found through their gamma-ray emission, 
long-wavelength surveys will soon provide a more powerful tool to pinpoint the nearest 
engine-driven supernovae. 
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Figure 1: Radio observations of the nearest massive star explosions. The 8.46 
GHz radio emission from SN 2009bb (red) is more luminous than any of the other 142 local 
{d < 200 Mpc) SNe Ibc observed (Ref.l2 and references within) to date on a comparable 
timescale {At < 100 days), and is consistent ^^'^^'^^ with the radio afterglow luminosities 
of the nearest GRBs discovered through their gamma-ray signal within a similar volume 
(black). Local SNe Ibc with well-studied radio emission (grey) exhibit lower luminosities 
and peak at later times, indicating smaller sizes and lower mean expansion velocities. The 
radio emission from most local SNe Ibc is below our current detection threshold; we include 
them here as upper limits {3a; grey triangles). 

Figure 2: Synchrotron self-absorption model fits to the SN2009bb radio spec- 
tra. The radio emission from GRBs and SNe Ibc is suppressed^^'^^ at low frequencies by 
SSA, defining the spectral peak frequency. The spectral shape below and above the peak 
is characterized as F,^ oc u^^"^ and F^, oc respectively, where p is the power- law 

index of the relativistic electron energy distribution above a minimum Lorentz factor, 7^. 
Our multi-frequency radio observations of SN 2009bb taken with the VLA and GMRT (see 
Suppl. Info.) on Apr 8, May 10, Jun 6-10, and Aug 8-11 UT {At w 20, 52, 81, and 145 
days) are well described by a standard SSA spectrum, with Up ^ 6, 3, 1, and 0.8 GHz, 
and peak flux densities of F,y^p f« 19, 15, 13, and 11 mJy, respectively. As the blastwave 
expands the shocked material becomes optically thin, causing Up to cascade'^^ to lower fre- 
quencies with time. The optically thin spectral index is constrained to roughly F^, ^ 
which implies p ~ 3, in line^^ with other radio SNe Ibc. Error bars are la. 

Figure 3: Radio properties of the neeirest massive steir explosions directly 
reveal the blastwave velocities. We compare the peak radio luminosities for SNe Ibc 
(red) and nearby GRBs {z < 0.1; blue squares) as observed at the spectral peak frequency. 
Up, and at time At. These observed properties are tightly related^^ to the blastwave radius. 
The average velocities are reasonably estimated as R/At (dashed grey fines). For SNe Ibc 
we infer typical velocities of R/At ^ 0.1c, while SN2009bb (star) and the nearest GRBs 
show R/At «i c. Error bars are la. 
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Figure 4: Blastwave velocity and energy for massive star explosions. We 

compare the blastwave velocities and inferred energies for the well-studied radio SNe Ibc 
and nearby GRBs included in Figure 3 with classic GRBs. Assuming equipartition of 
energy between electrons (eg) and magnetic fields (e^) as eg = = 0.33 which ac- 
counts for equal energy in shocked protons, the total internal energy of the radio emit- 
ting source is £^ = E mm /^R where the minimum energy, -Emmj is derived"^^ from the en- 
ergy density in magnetic fields, E^i^ ^ 8.3 x 10^6(//0.5)(S/G)2(i?/10^6 cm)^ erg, here 
B ^ 0.43(ee/eB) -4/19 (//0.5) -4/19 (L,,p/1028 erg s-^ Rz-^y^/^^Up/B GHz) G and / = 0.5 
is the fraction of the spherical volume occupied by the radio-emitting region. For radio SNe 
Ibc (red; Ref. 12 and references therein), the average velocities are reasonably estimated as 
/3 — R/cAt since the bulk flow is^^ expanding freely. The radio properties of these objects 
imply typical values of /3 ^ 0.1 and E ^ 10"^"^ erg. Adopting the same framework for 
SN2009bb (yellow star) we find R/At ^ 0.9c and E ^ 10^^ erg and note that the inclusion 
of relativistic effects (e.g., photon arrival time) would not change these results significantly 
as shown by Ref. 13 and 15 for GRB 980425/SN 1998bw. The nearest GRBs (blue squares) 
tend to be trans-relativistic and some decelerate rapidly, so we adopt the inferred^^'i^'i^'^^ 
blastwave velocities and energies from the earhest available radio data, at At < 5 days. 
Likewise, the inferred velocity of SN 2009bb at At ~ 20 days, is a strict lower limit since 
there is evidence that the blastwave decelerated early on (see Suppl. Info.). Finally, classic 
GRBs (blue circles) in the decelerating Blandford-McKee phase^^ have F oc t~^/^ so we 
conservatively estimate their blastwave velocity at At = 1 day according to F oc t~^/^ 
and adopt the beaming corrected blastwave energies from afterglow modeling (Ref. 6 and 
references therein). Error bars are la. 



A Relativistic SN 11 



T 1 — I — I I I I I I 1 1 — I — I I I I I I 1 — TT — I — I I I I I I \ : 

® GRB 980425 = 




1 10 102 

Time Since Explosion (days) 



Figure 1. 



12 Soderberg et al. 




Figure 2. 



A Relativistic SN 13 



io3or 




10 100 
(At/1 day)(i/p/5 GHz) 



Figure 3. 



14 Soderberg et al. 



1052 - 



>> 

1048 



© 



1-1 1 I r 

Gamma-ray Bursts 



SN 2009bb 









+ 



o 







IQ46\- Supernovae Ibc 



0.1 1 10 

Blastwave Velocity (Vp) 



Figure 4. 



Suppl. Info. 15 

SUPPLEMENTARY INFORMATION 
A Relativistic Supernova Powered by a Central Engine but 
Lacking a Detected 7— ray Burst 

1 Explosion Date 

The Chilean Automatic Supernova Search Program (CHASE; Ref. 31) and the Carnegie 
Supernova Project (CSP; Ref. 32) initiated a dedicated optical and near-IR follow-up 
campaign for SN2009bb beginning immediately after discovery. Using a combination of 
pre- and post-explosion data collected by CHASE and CSP a well-sampled optical light- 
curve has been constructed (G. P., M. D. S. & M. A. H. in prep). Fitting a expanding 
fireball model, L^pt oc t^, to the rising component of the optical light-curve implies an 
explosion date of to ~ Mar 18.4 UT ± 0.6 days. While this fit enables a first estimate of the 
explosion date, it assumes a constant velocity. Realistically some deceleration is expected, 
especially for fast expanding SNe hke SN 2009bb and SN 1998bw so it is reasonable that the 
explosion time should be somewhat later that the fireball fit. Along this hne, we match the 
early evolution to that observed^^ for SN 1998bw and find a slightly later explosion date, 
to ~ Mar 20 UT. We adopt an intermediate value of to = Mar 19±1 UT throughout the 
paper. Finally, we note that the inferred early discovery is consistent with the observed 
epoch of maximum light, Apr 1 UT, occurring two weeks after our estimated explosion 
date and in line with the observed rise times for other SNe Ibc (G. P., M. D. S. & M. A. H. 
in prep) . 

2 Radio Observations 

2.1 Very Large Array observations 

We observed SN 2009bb with the Very Large Array (VLA) on many epochs spanning Apr 

5.2 to Jun 6.0 UT. Data were collected at four frequencies, z/obs = 1.43,4.86,8.46 and 
22.5 GHz. All VLA observations were obtained in the standard continuum mode with 
2 X 50 MHz bands. At 22.5 GHz we included referenced pointing scans to correct for the 
systematic 10-20 arcsec pointing errors of the VLA antennas. We used 3C286 (J1331-I-305) 
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for flux calibration, while phase referencing was performed against calibrator J1036-3744. 
Data were reduced using standard packages within the Astronomical Image Processing 
System (AIPS). Flux density measurements were obtained by fitting a Gaussian model to 
the radio emission. For observations at 1.43 GHz we restricted the UV-range to reduce 
contamination from diffuse host galaxy emission from the strongly star-formation region 
in which the SN resides (Figure 5). The resulting flux density measurements are listed in 
Table 1 and shown in Figure 6. 

2.2 Giant Meterwave Radio Telescope Observations 

We observed the SN 2009bb with the Giant Meterwave Radio Telescope (GMRT) on several 
epochs spanning May 31 and Aug 10 UT at effective frequencies of 1287, 617, and 332 MHz. 
We used a 2 x 16 MHz bandwidth each divided into 128 frequency channels. The duration 
of the observations was two hours each, except for the 332 MHz observation which was 4 hrs 
in duration. We used 3C147 and 3C286 for flux and bandpass calibration and J1057-245, 
J1107-448, and J1018-317 to monitor the phase at 1287, 617, and 332 MHz, respectively. 

The GMRT data were reduced and analyzed with AIPS. To remove contamination from 
diffuse host galaxy emission, we removed the short baselines 0-5 kA and 0-3 kA at 1287 and 
617 MHz, respectively. At 332 MHz the host galaxy emission dominates that of the SN at 
the SN position. To estimate the flux density of the diffuse emission in the beam of the SN, 
we first construct a host galaxy map at 617 MHz and measure an integrated fiux density 
of Feu MHz = 113.5 ± 5.7 mJy after removal of the SN emission. We convolve this map 
down to the beam size of our 332 MHz map, using the AIPS task CONVL. The extended 
emission at the SN site is then scaled to 332 MHz using the spectral index {F^, oc of 
the host galaxy emission as measured between 1287 and 617 MHz. The 332 MHz map is 
then re-gridded to match the convolved 617 MHz map using the AIPS task HGEOM. The 
SN flux at the 332 MHz is then determined as the excess observed flux in the synthesized 
beam over and above the estimated extended emission at the SN position. Our resulting 
flux densities for all GMRT epochs are hsted in Table 2 and we note that the error in the 
SN flux density at 332 MHz is dominated by the uncertainty in determining the extended 
emission at 617 MHz. 
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3 Radio Temporal Evolution 

The temporal and spectral evolution of the synchrotron emission is determined by the 
dynamics, the density profile of the cjccta, and the properties of the circumstellar medium. 
Young, non-relativistic SNe expand freely since their ejecta are largely undecelerated while 
GRBs expand^^ according to the Blandford-McKee solution for relativistic, decelerated 
blastwaves. Both SNe and GRBs ultimately enter a phase of non-relativistic adiabatic 
expansion, the Sedov phase, when the ejecta have fully decelerated and the transverse 
velocity falls well below r/3 < 1. Trans-relativistic blastwaves, for which r/3 ~ a few, e.g., 
SN1998bw/GRB 980425, bridge these three hydrodynamic regimes and the flux density 
evolution depends on the structure of the ejecta. 

The early [At < 40 days) radio emission of SN 2009bb decays at both optically-thick 
and -thin frequencies (Figure 6), in contrast to the predicted^^ evolution for SNe Ibc in 
free expansion. This indicates that another process is responsible for powering the ra- 
dio emission at early time: a central engine. In this context, the early evolution may be 
understood if the geometry of the relativistic outflow was originally aspherical. In this 
scenario, the relativistic ejecta are detached from the freely expanding SN outflow and 
decelerated^ on a timescale, idee ~ 1 {E/IO"^^ erg)(M/10-^ Mq yr"^)"^ day, when r/3 < 1. 
Here, M — Anpy^r'^ is the mass loss rate of the progenitor star, p is the circumstellar 
density and is the velocity of the progenitor wind. Based on the SSA model fits to the 
radio spectra, we estimate^ the mass loss rate of the SN 2009bb progenitor to be M w 6.4 x 
10-^(ee/eB)-«/^9(eB/0.33)-i(//0.5)-«/i9(F,,p/mJy)-Vi9(^^/5 GRzf {At/10 days)^ yr"! 
or M 2.0 X 10-^ Mq yr"^ for / = 0.5 and eg = = 0.33 at At ^ 20 days. Here we have 
assumed a stellar wind density profile, p oc r~^, and a steady progenitor wind velocity of 
Vw ~ 10^ km s~^ consistent with the observed^^ properties for local Wolf-Rayet stars. For 
SN 2009bb, we estimate that the outflow decelerates by tdec ~ 6 days and note that this 
estimate is independent of the partition fractions. Therefore the mean expansion velocity 
of R/ At ~ 0.9c inferred at At 20 days is a strict lower limit on the initial velocity of 
the outflow. Thereafter the ejecta spread^^ sideways, stalling the radial expansion as the 
outflow becomes spherical and causing the peak radio emission to decay. The blastwave 
spherizes^^ on a timescale, tsph ~ 6 {E/IO'^^ erg)(M/10-^ M© yr-^)-^ days. By At ^ 34 
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days, we estimate that the SN2009bb blastwave was spherical and the radio evolution 
thereafter is well described by a standard spherically symmetric model. 

4 X-ray Observations 

4.1 Chandra data 

We observed SN2009bb with Chandra ACIS-S beginning at Apr 18.6 UT {At ^ 31 days) 
for 10 ksec under our approved Cycle 10 program to study the X-ray properties of radio- 
emitting Type Ibc supernovae. We detect an X-ray source coincident with the optical and 
radio counterparts at position q;(J2000)=10''31'^33.86' and 5(J2000)=-39°5r29.4" (±0.4 
arcsec in each coordinate) . The data were reduced and the source counts were extracted in 
the standard manner using the Chandra threads. Within a 4.9 arcsec aperture we extract 
24 counts. After correcting for an inferred foreground extinction of E{B — y) fa 0.5 (see 
Ref. 40) corresponding^^ to a column density, uh — 2.92 x 10^^ cm~^, we estimate an 
(absorbed) unabsorbed flux of {Fx = 1.5 ± 0.3 x 10^^ erg cm'^ s"^) Fx = 2.3 ± 0.5 x 
10~^^ erg cm~^ s~^ (0.3-10 keV) assuming a power-law spectrum with photon index, F = 2. 
At the distance of NGC 3278, the unabsorbed flux corresponds to an X-ray luminosity of 
Lx — 4:A ± 0.9 X 10^^ erg s~^. The X-ray emission from SN 2009bb is therefore a factor of 
a few less luminous than GRB 980425 observed on a comparable timescale and at the high 
end of the X-ray luminosities observed^^ for other local SNe Ibc. 

4.2 Swift X-ray Telescope Data 

We also observed SN2009bb with the Swift X-ray Telescope (XRT) on several epochs 
spanning Mar 24.1 to Apr 18.7 UT. Data were retrieved from the HEASARC archive and 
reduction performed with the xrtpipeline script packaged within the HEAsoft software. 
We used the default grade selections and screening parameters to extract measurements 
for the 0.3-10 keV energy range. 

We extracted the counts within a 47 arcsec (90% PSF containment) aperture centered 
on the optical SN position in each epoch. Diffuse emission from the host galaxy is clearly 
detected within the extraction region. We bootstrap the flux calibration to the Chandra 
data, temporally coincident with our final XRT epoch. Adopting the same 47 arcsec 
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aperture for the Chandra data, we extract 233 counts after subtracting off the background 
count rate. Fitting the same absorbed power-law model described above, we find an 
unabsorbed flux of Fx = 3.0 ± 0.8 x 10^^^ erg cm'^ s'^. Therefore, within the XRT 
PSF, the host galaxy emission dominates that of the SN by a factor of ~ 10. Since the 
fractional uncertainties of the XRT fluxes preclude our ability to extract estimates of the 
SN flux, we adopt the measurements as upper limits (see Table 3). Finally, we note that 
a comparison of the first XRT upper hmit and the measured SN flux from the Chandra 
observation constrain the temporal evolution of the SN X-ray flux to be no steeper than 
Fx oc t-\ 

The flux ratio between radio and X-ray band can reveal the nature of the X-ray emission 
which, in the case of SNe Ibc, may be dominated by (i) thermal emission associated with 
the shock-heated ejecta, (ii) non-thermal synchrotron emission from the shocked circum- 
stellar material, or (iii) inverse Compton up-scattering of optical photons by radio-emitting 
electrons (see Ref. 37 for a review). In the case of SN2009bb, at the epoch of the Chan- 
dra observation we measure the spectral index between the radio {ur — 8.46 GHz) and 
X-ray band (i/^ = 2 x 10^^ Hz), to be F^ oc i/-o-82±o.o2 rpj^-g -g somewhat flatter than 
that expected by an extrapolation of the radio synchrotron spectrum as F^, oc Given 
that the optical SN emission peaks on a similar timescale, it is feasible that the X-ray 
emission is dominated by inverse Compton emission. Alternatively, if the shock is cosmic 
ray dominated the synchrotron spectrum flattens^^ at higher frequencies. The nature of 
the X-ray emission will be discussed in a subsequent paper. 

5 Searches for a Coincident High Energy Signal 

Massive star explosions can give rise to high energy emission at the moment of explosion 
via two primary channels: (i) strong X-ray and/or gamma-ray emission powered by ultra- 
relativistic ejecta as inferred for GRBs, and (ii) weak X-ray emission produced by shock 
breakout in ordinary non-relativistic supernova explosions (e.g., SN2008D). In order to 
search for a statistically significant association between an optical SN and a high energy 
transient, an accurate estimate of the explosion time is required. SN2009bb is unique 
among local optically-discovered SNe in that the optical data set is sufficient to enable a 
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reasonably tight constraint on the explosion date (§1). Along this line, we searched for 
high energy counterparts between the dates Mar 18-20 UT and coincident with the SN 
position. 

Among the automatic triggered public bursts announcements made through the 
Gamma-Ray Burst Coordination Network (GCN) from Swift, Fermi, Integral, Agile, 
Suzaku and the Interplanetary Network (IPN), we note just one marginally coincident 
burst. The weak Fermi Gamma-ray Burst Monitor (GBM) burst GRB 090320C was de- 
tected on Mar 20.045 UT at q;(J2000) = 07*^13"^, 5(J2000) = -43°18' {la positional uncer- 
tainty of 15 deg; Ref. 42) and lies 3.1a away from the optical and radio SN 2009bb positions. 
With a duration of Tqq 4 sec, the flux of the burst is = 3.9 x 10~^ erg cm~^ s~^ (8- 
1000 keV). Given the large positional uncertainty due to the weakness of the burst, we 
undertook a Monte Carlo experiment to calculate the probability that GRB 090320C was 
associated with SN2009bb. We used 1000 iterations to test the frequency that a GRB 
with the properties observed for GRB 090320C could be recovered from a burst within la 
of the location of SN2009bb given the known detector response. From this experiment, 
we exclude the association of GRB 090320C and SN 2009bb at the 99% confidence level. 
Separately, we note that the SN position was in the field-of-view of the GBM for 62 percent 
of our bracketed explosion date estimate. 

We next searched for sub-threshold bursts in the Swift Burst Alert Telescope (BAT) 
data, including those data from the BAT Slew Survey (BATSS; Ref. 43) for the same 
range of dates. The SN position was within the 50 deg field-of-view of the BAT for just 7 
percent of this time. With an accuracy of ~ 2 arcmin, no sub-threshold (Fy > 5.5a above 
background) triggers were detected at a position coincident with SN2009bb. Moreover, 
in the BATSS data set we find no statistically significant {Fj > 3.5a) bursts during the 
spacecraft slews over this timescale. We also note that GBM burst 090320C was not in 
the field-of-view of the BAT at the time of the Fermi trigger. 

By combining the monitoring times of BAT and GBM and accounting for the 9772 
seconds when the SN position was monitored by both satellites simultaneously, we estimate 
a combined monitoring time of 65 percent. Furthermore, both GBM and BAT are included 
within the all-sky Interplanetary Network (IPN; ) of gamma-ray satellites which offers 100% 
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temporal coverage to SN 2009bb over the bracketed explosion date estimate. Based on the 
detection efficiency of IPN as reported in Ref. 44, we place an upper limit on the detectable 
gamma-ray fluence from SN2009bb of < 5 x 10"^ erg cm-^ (band, 25-150 keV). At 
a distance of d «i 40 Mpc, this corresponds to an isotropic-equivalent energy release of 

< 10''^ erg which is a factor of two higher than that observed^^ from GRB 980425. 

The lack of a detected gamma-ray counterpart can be attributed to several factors. 
For example, if the GRB jet is not directed towards our line-of-sight, relativistic beaming 
inhibits the detection of gamma-rays. Most importantly, as seen in the case of the nearest 
GRBs, the gamma-rays can be exceedingly weak and the isotropic gamma-ray energy yield 
dwarfed by that of the blastwave kinetic energy (e.g., GRB 980425; Ref. 45). Along this 
line we note that the rate of such intrinsically weak GRBs at 2; < 0.1 is 3-10 times higher 
than that of classic GRBs discovered preferentially at larger redshifts. 

6 Predictions for the Rate of Engine-driven SNe from Long- Wavelength Sur- 
veys 

Local {d < 200 Mpc) SNe Ibc are currently discovered* at a rate of ~ 25 yr~^ by the 
ongoing efforts of amateur astronomers and automated searches (e.g.. Lick Observatory 
Supernova Search, CHASE) that monitor weekly a sample of nearby galaxies searching 
for new or transient objects. With the advent of wide-field optical surveys, specifically 
The Palomar Transient Factory (PTF; Ref. 47) and Pan-STARRS,^^ this sample will be 
supplemented by similarly nearby and young SNe Ibc at a rate of ~ 50 yr~^ (PTF) and 
~ 25 yr~^ (Pan-STARRS) over the next 3-5 yrs given the volumetric rate^^'^° of SNe Ibc 
of 1.7 X 10^ Gpc-V~^- 

Through prompt radio observations of the 100 SNe yr~^ discovered by these combined 
efforts, engine-driven events like SN 2009bb can be clearly identified by their luminous and 
low-frequency radio spectra. The current sensitivity of the Very Large Array {3a ^ 0.1 
mJy in 15 min) implies that the radio luminosity hght-curve of SN 2009bb extending to 
~ 100 days could be detected with S/N > 5 out to a distance of 130 Mpc. The increase 
in sensitivity of the Expanded- VLA (expected 2010; Ref. 51) will enable such events to 
*http://www.cfa.harvard.edu/iau/lists/Supernovae.html 
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be detected significantly further, to 400 Mpc, well matched to the volumes probed by the 
optical surveys, and requiring just 25 hrs of observing time per year. Given our estimated 
fraction of engine-driven events, 0.7t.o^%, we anticipate that ~ 1 new relativistic SN will 
be both discovered and identified each year. This is a factor of three higher than the 
current discovery rate^^ of nearby GRBs with gamma-ray satellites, demonstrating that 
long-wavelength surveys will ultimately become powerful tools in the discovery of such 
events. In fact, over the next decade, the "deep-wide-fast" optical survey planned for the 
Large Synoptic Survey Telescope (LSST; Ref. 52) will reveal ~ 500 SNe Ibc each year 
within 400 Mpc, of which up to ~ 10 will be engine-driven. 
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Date 




1^,4. OD t^rlz 


^ /^,o.40 Lrrlz 


^ l/,ZZ.O ijxlZ 


VLA 




f niJvl 


fniJv) 


fniJv) 


f niJv) 


Config. 


Apr 5.2 






24 681 ± 066 




B 


Apr 8.2 


1 308 ± 099 


12 863 ± 161 


17 568 ± 088 


5 204 ± 191 


B 


Apr 13.2 






16 349 ± 107 




B 


Apr 15.1 




9.629 ± 0.181 


13.812 ± 0.114 


5 170 ± 394 


B 


Apr 23.2 


1.141 ± 0.159 


7.601 ± 0.226 


8.881 ± 0.121 




B 


Apr 29.1 


2 963 ± 348 


7.632 ± 0.202 


7 714 ± 095 


1 796 ± 240 


B 


Mav 3 1 


3 406 ± 122 


9 680 ± 165 


8 482 ± 098 




B 


Mav 10 1 


4 939 ± 132 


12 215 ± 196 

J. ^ • ^ J. C_/ 1 \JmJ.lJ\J 


6 824 ± 102 




B 


Mav 13 


5 836 ± 136 


9 727 ± 508 


6 327 ± 151 




B 


Mav 20 1 


6 903 ± 1 80 


8 542 ± 1 72 


3 294 ±0118 




B 


May 27.0 


6.702 ± 0.351 


8.568 ± 0.080 


4.204 ± 0.060 




CnB 


Jun 6.0 


9.465 ±0.186 


5.877 ±0.110 


3.203 ±0.074 




CnB 


Jun 17.0 


12.133 ± 0.249 


4.220 ± 0.120 


2.392 ± 0.082 




CnB 


Jun 26.0 


12.189 ± 0.522 


3.556 ± 0.350 


1.903 ±0.548 




C 


Jul 18.9 


11.230 ±0.356 


2.510 ±0.161 


1.032 ±0.104 




C 


Aug 11.8 


8.687 ± 0.341 


2.477 ±0.155 


1.084 ±0.091 




C 



Table 1. Very Large Array observations of SN2009bb. 



Date 
(UT) 


Frequency 
(MHz) 


Flux Density 
(mJy) 


Synthesized Beam Size 
(arcsec) 


May 31.5 


1287 


4.6 ±0.3 


5x2 


Jun 10.6 


617 


3.0 ±0.3 


10x4 


Aug 8.3 


617 


7.9 ±0.3 


10x4 


Aug 10.4 


332 


3.3 ±0.7 


13x7 



Table 2. GMRT observations of SN2009bb. 



Date 


Exposure 


Flux 




Mission 


(UT) 


(sec) 


(erg cm~^ s' 






Mar 24.1 


3964 


< 1.3 X 10- 


-13 


Swift/XRT 


Apr 6.2 


9829 


<1.7x 10- 


13 


Swift/XRT 


Apr 10.5 


6720 


< 2.5 X 10- 


-13 


Swift/XRT 


Apr 18.7 


12752 


< 3.0 X 10- 


-13 


Swift/XRT 


Apr 18.6 


9914 


2.3 ±0.5 X 10-1^ 


Chandra 



Table 3. X-ray observations of SN2009bb. All measurements are 0.3-10 keV and unabsorbed. 
Upper limits are dominated by diffuse host galaxy emission. 
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Figure 5. A composite image of the host galaxy of SN2009bb, NGC3278, constructed from 
2MASS near-IR (red), Swift/VYOT optical (V-band, green) and UV (UVW2-filter, blue) images. 
Contours defining the structure of the diffuse radio (1.43 GHz) emission as observed with the 
VLA in B-configuration are over-plotted and trace the regions of star-formation in the galaxy. 
The SN is clearly detected in the radio and optical images (labeled) and lies within a strongly 
star-forming region of the galaxy. 
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Figure 6. Radio light-curves for SN2009bb at multiple frequencies. 



